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8. Evidence from Spectra of Bright Fireballs 



ZPENEK CePLECHA 

AstvoYioniical Institute of the Czechoslovak Academy of Sciences 

Ondfej ov, Czechoslovakia 


Spectral data with dispersions from 11 to 94 A/mm on 4 fireballs of actual brightness 
of —4 to —12 magnitude and with velocities of about SO km/s at 70 to 80 km heights 
are used for studies of meteor radiation problems. Previously published analyses need 
revision for two main reasons: (a) the absolute values of oscillator strengths of Fe I lines 
from laboratory data were recently recognized to be 1 order of magnitude lower, (b) the 
luminous efficiency factor r of Fe I is now much better known from several different 
experiments. The radiation of fireballs is found to be strongly affected by self-absorp- 
tion. But if the emission curve of growth is used for correction of the self -absorption of 
Fe I lines, a great discrepancy between spectral data and efficiency data for total Fe I 
light is found. If one assumes that the self -absorption is superposed on another effect, 
a decrease of the dimensions of the radiating volume with increasing lower potential Ei, 
the spectral data on Fe I lines will be in agreement with the luminous efficiency of total 
Fe I meteor radiation. Formulas for emission curve of growth and Boltzmann distribu- 
tion including this effect are derived. This effect is important for fireballs brighter than 
about -1 or -2 magnitude, while self-absorption seems to be important even for fainter 
meteors. The optically thin radiation of all Fe I Urns might be expected for meteors 
fainter than +5 magnitude. Excitation temperature of 5500° K and relaxation time of 
0.02 s were found as typical values for the Fe I radiation of fireballs studied . The hg t 
of fireballs is emitted during a relatively long relaxation time, which is many orders of 
magnitude longer than the time necessary for spontaneous radiation of excited Fe 
atoms. The dimensions of the radiating volume of Fe I gas for lines with Ei =0 were 
found tobe 0.3X9 m at 0 absolute magnitude and 2X60 m at —10 absolute magnitude. 
It was not possible to determine any realistic abundances of other elements due to sma 
numbers of lines for an analysis independent of Fe I, while the Fe I curve of grout 
cannot be need for other elements, because the radiation originates mainly from the 
effective surface of the radiating volume. *4 general formula for meteor radiation is also 
derived and compared with the conventional luminosity equation. 


T here is enough light available to obtain 
spectra of bright fireballs with dispersions 
well below 100 A/mm, On the other hand, there 
arc not so many bright fireballs which is, in part, 
the reason that published spectral data on line 
intensities arc rather scarce. Preparing this paper 


I used, almost exclusively, my own observational 
material resulting from ten years of systematic 
operation of three grating spectrographs employed 
in the registration of meteors at the Ondrejov 
Observatory in Czechoslovakia. On an average, 
about one year was necessary to get a good spec- 
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trum of a fireball (over 100 spectral lines in the 
visible region) using these cameras. The direct 
double-station cameras (30 in number) provided 
us with all the necessary geometrical and velocity 
data on these fireballs. The reduction of four 
spectra is now completed; three of them have been 
published already (Ceplecha and Raj chi, 1963; 
Ceplecha, 1964, 1965, 1966) and one spectrum, 
S 1132, containing 990 spectral lines, is in print 
(Ceplecha, 1971). 

The summary of the spectra I am using in this 
study is in table 1. All the objects were of approxi- 
mately the same velocity, about 30 km/s. Thus 
the observational data presented here are to be 
referred to this particular velocity, to heights of 
70 to 80 km and to an actual brightness interval 
of —4 to — 12 absolute magnitude (the maximum 
absolute magnitudes are in the interval from —9 


to - 12 mag). I also used qualitative data on one 
more spectrum, belonging to a fireball of approxi- 
mately — 18 magnitude. The reduction of this 
spectrum is not yet finished and only preliminary 
data are available. The spectrum belongs to 
a body with a significant terminal mass of about 
600 g, but no meteorite was found. 


NEED FOR REVISION OF PREVIOUS 
STUDIES 

Some studies of spectra in table 1, which I 
published previously, need a revision from two 
main aspects: 

(1) The absolute values of the oscillator 
strengths of Fe I lines from laboratory data were 
recently recognized to be one order of magnitude 


Table 1—List of Meteor Spectra Used for A nalysis in This Paper 


Spectrum No. 


Meteor No. 
Date 

Time (GMT) 

a M rw n 


Appearance on 
spectral plate 


Spectral order 
Dispersion (A/mm) 


No. of spectral lines 
«« (deg) 

(deg) 
cos Zr 
v*> (km/s) 
a (AU) 
e 

<7(AU) 

a? (deg) 1 

(deg) h 950.0 
* (deg) I 


> Mn: 


S 6 

S 526 

S 912 

130960 

32281 

36221 

Sept. 13, 1960 

Nov. 13, 1961 

Sept. 26, 1962 

18 h 59" 

21 h 49 ra 46 a 

23 h 18 m 16“ 

-30.5 

-9.6 

-10.3 

flare of 0.025 s 

trail of 0.6 s with 

trail of 0.4 s with 


many flares 

steep but smooth 
light curve 

2 to 7 

1+2 

1 

38 to 11 

94+47 

58 

166 

161 

189 


58.2 

22,2 

~0 

16.5 

7.1 

«0.2 

0.785 

0.718 

~30 

28.86 

31.92 


2.30 

1.63 


0.814 

0.849 

Southern Taurid 

0.427 

0.247 


105.5 

129.9 


51.1 

3.2 


5.0 

4.9 


S 1132 


(S 1135) 


i 


38421 

Apr. 21, 1963 
l h 13 m 48® 
-12.4 


flare of 
0.022 s 


1 

57 


trail of 0.4 s 
with 
terminal 
flare 


pnsm 


990 

264.6 

56.8 

0.965 

32.57 

5.51 

0.819 

0.997 

190.8 

30.0 

49.0 


is the total absolute maximum magnitude from panchromatic plates, while magnitudes in table 2 are the ac- 

mnfrmttirlAa n ^ C TT^ T I ! _ i ! 1 4 


tual absolute magnitudes at each point of Fe I radiation only. 



EVIDENCE FROM SPECTRA OF BRIGHT FIREBALLS 


91 


lower than previously accepted (Grasdalen et ah, 
1969; Garz and Kock, 1969) . 

(2) The luminous efficiency factor r of Fe I is 
now much better known from several different 
experiments (Ayers et al. 1970, Friichtenicht 
et al. 1968). 

The absolute abundances of Fe I derived from 
spectral data on fireballs were previously recog- 
nized by me to be about 1 or 2 orders of magnitude 
greater than that derived from total Fel light, 
using a luminous efficiency factor (Ccplecha, 1965, 
1968). The “spectroscopic” mass of meteors 
studied is increased by an additional order of 
magnitude over the “photometric” mass, if we use 
the new Fe I oscillator-strength data. The main 
purpose of this paper is a comparison of the spec- 
tral data on Fc I absolute abundances with 
corresponding data derived from the total Fe I 
light curve calibrated by the luminous efficiency 
factor. 

THE LUMINOSITY EQUATION 

Before dealing with the main problem of this 
paper, I will derive a new 7 general form of lumi- 
nosity equation, which wall help to relate the 
conventional luminosity equation to the gas mass 
in the radiating volume. 

Let us denote by tyia the total gas mass that has 
entered the radiating volume since the meteor 
beginning. Let us denote by mz the total gas mass 
that has left the radiating volume since the 
meteor beginning. Then, at each instant, the total 
gas mass vig inside the radiating volume is 

mg — Wa — ( 1 ) 

Assuming that the excitation temperature and 
self-absorption are roughly constant during the 
meteor flight, which is applicable to radiation of 
one clement in one state, w 7 e have for the intensity 
I of the radiation : 

I = km G (2) 

where the proportionality factor k is a luminous 
efficiency of the radiating gas measured in 
ergs s _I g -1 . Even if k is some function of time, it 
changes substantially less than tyig and the propor- 
tionality, equation (2), may be considered a good 
approximation to the reality. This approximation 


is no worse than the conventional luminosity 
equation (4). 

Substituting equation (1) in equation (2) we 
have 

I = k (ni A — rriz) (3) 

This is a meteor luminosity equation more general 
than the conventional equation 



Our task now 7 is to relate the general equation 
(3) to the special equation (4). The additional 
assumptions to be applied to equation (3) can be 
derived from the following considerations: We 
differentiate (3) 

dl = k (dwiA — dmz) (^) 

Wc can assume that dm a and dmz are the same 
functions of time with different phase only 

dmz(t)=dm A (t-t R ) ( 6 ) 

where t R is constant: it is the time spent by an 
average atom inside the volume of radiating gas, 
and I will refer to it as the relaxation time in this 
paper. Now 7 if equation (6) is substituted in 
equation (5) and Taylor's series with the first 
derivative used, we have 

dl = Mr d (dm a / di) (7) 

If we replace m A by the ablated meteor mass 

m A — m x —m (8) 

and integrate, we have 

I=—kt R {dm/dt) +k R (9) 

But there is practically no light at the beginning 
of the ablation of the meteor and k R = 0. Thus 
equation (9) is equivalent to equation (4 ) , if 


kin** 


rf_ 

2 


( 10 ) 


Variation in v and r during the trajectory of one 
meteor is neglected. We can summarize: 

(1) The conventional luminous equation (4) is 
only a very 7 special case of the more general 

equation (3). , 

(2) The conventional luminous equation (4) 
represents an average smooth change of light. 
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(3) Equation (3) can easily explain the big 
fluctuations of the observed light curves (see 
fig. 1) with dm/dt being a smooth function of 
time. If dm z < dm a , I increases; if dm z >dm Ay I 
decreases. Thus dm a , dm z may be smooth func- 
tions of time, while the light intensity is flaring. 
Figure 1 shows this clearly even for Ir— const.; 
but without this limiting assumption, the general 
form of equation (3) is even more sensitive for 
producing big changes in light emission from small 
smooth changes of dm/dt . 

SELF-ABSORPTION IN SPECTRA OF 
FIREBALLS 

I have already published studies of three spectra 
of meteors in table 1 using the method of emission- 
curve-of-growth, which takes into account the 
optical thickness of individual spectral lines 
(Ceplecha 1964, 1965, 1967). After publication of 



Figure 1 . — Theoretical light curve computed from arbi- 
trarily chosen but smooth values of dm a, dmz. This 
shows the flaring effect. Formula (3) was used with 
simplifying assumption (6): log l R = -2.7, t» = 30 
km/s, log £ = 3.42. This is just an example represent- 
ing formula (3) with limitation (6). General form of 
formula (3) is even more able to describe the flaring 
effect. 


this method, there were some doubts as to whether 
the method is applicable to fireball spectra or not. 
The importance of self-absorption and damping 
processes for fireball radiation may be well demon- 
strated using figure 2. Even though figure 2 is 
presented here only for spectrum S 6, the same 
relative positions of individual multiplets result, 
using any of the other 3 spectra of table 1 at any 
of the measured points. The upper level potential 
E 2 is plotted along the x-axis and the total number 
of atoms of Fe I in all levels of the upper term, 
jV(term) , is plotted along the y-axis. The observed 



Figure 2. Numbers of Fe I atoms, iV(tcrm), occupying 
all levels of the upper term were computed for B 6 
spectrum of a flare on assumption that the radiation 
observed in corresponding lines was optically thin 
= log iV(term) is plotted against the upper po- 
tential Ei. A straight line corresponds to Boltzmann 
distribution of levels. One average value for each 
multiplet is given. The numbering of the points cor- 
responds to the notation of multiplets after Moore 
(1945). A self-absorption for multiplets with the same 
upper level term is evident: multiplet 4 against 15, 
5 against 16, 20 against 62, 21 against 39, 41 against 
168. The parameter Ir is the relaxation time, if the 
duration of the flare. The extrapolation of the Boltz- 
mann straight line to #2 = 0 must be equal to or lower 
than the point l R =l / . 
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radiation is assumed to be optically thin here, 
which means 


log N (term) 


ZoX 3 m e G 2 
8t rhehjf 


(ID 


Jo is optically thin radiation: here the observed I 
was used for figure 2 instead. G 2 is the total 
statistical weight of the upper term. Each point 
in figure 2 corresponds to an average value for the 
whole multiplet assuming that no self-absorption 
is present ( Iq=I ). Some of the multiplets with 
more observed lines point to a dependence of 
log N (term) on log gf as is apparent from figure 3. 
In figure 2, multiplets with the same upper level 
term yield greater log 2V (term), if their lower 
level is higher. Thus multiplet 4 (2?i = 0) shows 
an absorption of more than one order of mag- 
nitude relative to multiplet 15 (£i = 0.9 V), 
multiplet 5 more absorption than multiplet 16, 
multiplet 20 more than 62, multiplet 21 more than 
39 and multiplet 41 more absorption than multi- 
plet 168. This clearly means that we sec only a 
part of the light produced inside the radiating gas 
volume and emitted in multiplets 4, 5, 20, 21, 41. 
From other spectra than S 6, multiplets 2, 23, 


42, 43 and many others show clear evidence of 
self-absorption. But there is no reason for the 
assumption that multiplets 15, 16, 62, 39, 168, 
giving the greatest population observed for the 
same upper level, are optically thin. If we visualize 
a Boltzmann distribution of levels in figure 2, we 
have a straight line, the slope of which corre- 
sponds to the excitation temperature T: 

5040 

log N (term) = — 0 # 2 + const., where 0 — ^ 

( 12 ) 

The extrapolation of this straight line to E 2 = 0 
gives the total number of atoms in the ground 
state, which is almost the number of all atoms of 
Fe I for all reasonable temperatures. On the other 
hand, spectrum S 6 belongs to a short bright flare, 
the duration of which (0.025 s) combined writh 
the luminous efficiency of Fe I (Ayers et al., 1970, 
log T pa *= — 11.94) gives an upper limit for log N . 
This is equivalent to an assumption that the 
duration of the flare is identical with the relaxation 
time i R , If the relaxation time is less, the limit wall 
be lov'er and the discrepancy higher. From figure 



Figure 3.— Log jV(ierm), from individual lines of multiplets 4 and 15 (upper term z'FP), is 
plotted against log gf of the lines, assuming optically thin radiation as in figure 2. Self- 
absorption is evident. The dotted lines are only schematic, but indicate that lines with log 
gf —5 and less might be expected optically thin for z 5 D° of this spectrum. 
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2, it is also evident that multiplets giving the 
greatest log iV(term) must be affected by some 
further absorption effect. As a byproduct we can 
guess that the relaxation time t H of Fe I lines with 
Ei = 0 is of the order of 10 -2 or 10 -3 s, if reasonable 
excitation temperatures are used. There is no 
reasonable way to make t R less than 10 -4 s. This 
time is by orders of magnitude longer than the 
time necessary for spontaneous radiation of the 
corresponding excited states of the Fe I atoms. 
Some quasiequilibrium conditions may be reached 
during this relatively long time. There is no 
doubt, inspecting figure 2, that self-absorption is 
at least partly responsible for the observed effects. 

But if we use the emission curve of growth to 
correct for self-absorption, we get a relaxation 
time, bf, still greater than the duration of the flare. 
I tried to use a model with exponential decrease of 
temperature along the radiating volume to avoid 
this discrepancy, but it was not much help. This 
suggests that another additional effect is present. 

The observed radiation and deduced number of 
atoms are related to the whole volume of radiating 
gas. As we know that the wake phenomenon is 
strongly dependent on the excitation, we may 
expect the same dependence also for the “head” 
radiation. Then wc actually do see the radiation of 
different spectral lines belonging to different 
volumes of radiating gas with different relaxation 
times. The low potential lines perhaps radiate 
from a larger volume with longer relaxation time. 

I will present here a simple theory of this effect 
and derive corresponding changes of the emission 
curve of growth and of the Boltzmann distribu- 
tion. I will apply the theoretical results to the 
spectra of table 1 . 

EMISSION CURVE OF GROWTH WITH 
RESPECT TO CHANGES OF THE VOLUME 
OF RADIATING GAS WITH THE E 1 
POTENTIAL 

I will use the same simple model of radiating 
gas volume that I published previously (19G4), 
which is a cylindrical volume with great density 
on the surface, termed the "cylindrical shock- 
wave” model. I assume that the length t and 
radius R of the cylindrical volume are propor- 
tional, p being the proportionality factor. Then 

t=pR (13) 


The intensity of the wake phenomenon seems to 
depend on the lower level E x of the corresponding 
line. Thus I assumed a dependence of l and R on 
the lower level E x having an exponential form 

1{E X ) =/(0) exp l-2q Xl /{h B T)^ (14) 

where xi is the energy of the lower level (potential 
Ei) and q is an exponential factor describing the 
decrease of the length of the radiating volume with 
increasing lower potential. Here f(0) is the length 
of the radiating volume corresponding to the 
ground state. Following equation (18) of my 
previous paper (1964), we have 

log N(Ei) =log N(0) -5qQE 1 (15) 

for the total number of all atoms of the ground 
state inside a volume corresponding to the E x 
level. Here N (0) is the total number of atoms of 
the ground state inside the radiating volume 
belonging to the ground state. 

The Boltzmann distribution for this ease is then 
represented by 

log =~&(E 2 +5qE 1 ) +Iog — — +log iV(0) 
yj vYi e u 

(16) 

where 7 0 is the optically thin intensity radiated by 
the assumed changing volume in all directions. 
The observed intensities I arc to be converted to 
h by using the emission curve of growth given in 
our partial case by 

l°g ~ +4 qBE l 


+log 


c 

&rt; m [fi*(0)+fl(0);(0)] cos a 


=log F(X) 


(17) 

where 

log X = log gfa — GEi - qBEi 

. 1 3eW(0) 

+iog — nsi 

& mrn e cUv m ZR 2 (0) +R{ 0)7(0)] cos a ' * 

and the function F(X) describes the actual shape 
of the curve of growth. This is the same function 
as previously published [Ceplecha, 1964, equation 
(22)]. The damping part of the emission curve of 
growth is also well described by 
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lim log F ( X ) = 3^2 l°g X + const. (18a) 

in this case and generally. 

The procedure used for application of the above 
derived formulas to the observed spectral lines is 
exactly the same as published previously (Ceplecha, 
19G4). The results of the computations arc: the 
excitation temperature T, the number of all 
atoms of one clement in one state inside a radiating 
volume corresponding to Ei = O f and I n m propor- 
tional to the surface of this volume. A set of these 
values is obtained for each chosen value of q. If 
the ratio p of l to R is assumed, the volume of 
radiating gas for Ei = 0 is determined and the 
terminal density no(5o) at the end of the cylindrical 
volume as well as the average density n( 5) may be 
computed. The whole procedure of numerical 
computations was programmed for the Minsk-22 
digital computer. 

We have actually added one more parameter, q f 
to be determined from observation. But this 
parameter may be guessed from the spectia of 
flares of short duration. If dm/di is determined 
from equation (4), and compared with the results 
from the emission curve of growth, we can deter- 
mine the relaxation time 


mp 

dmjdi 


m G = ilN 


(19) 


Here m<? is determined from the curve of growth 
and dmjdi from the light curve. We have an 
additional condition for short duration flares 


t R <t f (20) 


where tj is the duration of the flare. 


RESULTS FOR FE I 

The actual solution of the problem was possible 
for Fe I lines almost exclusively. Any other 
element, except for Na I, Ca I and 0 I in spectrum 
S 1132, had not a sufficient number of spectral 
lines for such a solution. Only the lines without any 
blends were used for all the computations. The 
{//-values used were from Corliss and Bozman 
(19G2) and Corliss and Warner (1964), but 
corrected by —1.00 in the log to the newly 
established approximate absolute scale (Grasdalcn 
et ah, 1969; Garz and Kock, 1969). The excitation 
potentials were taken from Moore ( 1945) . 


I started with S6 and S 1132 and found a 
temperature of 5500° K as the lowest possible to 
satisfy the equation (20) condition in both 
spectra. With this temperature, t H of S 6 is equal 
to the duration of the flare of S 6, and t R of S 1132 
is about H nf the duration of the flare of S 1132. 
Any temperature higher than 5500° Tv may be a 
solution of the problem as well, but there is 
another rough estimation possible using the com-* 
puted density inside the radiating volume, the 
terminal Ro should not exceed too much the 
density of the surrounding air. If n 0 of Fe I is not 
greater than 10 times the surrounding air, the 
temperature should not exceed some 8000° K. 
If n 0 is of the order of the surrounding air, we 
arrive at about 5000° to 6000 Iv. Thus ooOO Iv 
was taken as a standard temperature and all the 
solutions for q were made with this temperature. 
A graphical example of one solution for S 1132 
is given in figures 4 and 5. The curve of growth, 
figure 4, and the Boltzmann distribution, figure 
5, are used simultaneously for the temperature 
determination. This means that the actual posi- 
tion of the theoretical curve of growth and the 
Boltzmann straight line are dependent on the 
positions of points in both figures 4 and 5. 

All points of the meteor spectra of table 1 were 
used for the solution of q for Fc I lines, using the 
above estimated temperature. The following 
values were determined and are listed in table 2; 
number of all atoms of Fe I inside the radiating 
volume; Y Um proportional to the surface of the 
radiating volume; dmjdi from the light curve of 
total Fe I radiation (log r pa n= — H*94 was used 
for I being in 0 magnitude units combined with 
e.g.s.; Ayers et al., 1970); the ratio of Fe I 
radiation to total radiation; the relaxation time 
t R \ rough guesses for the ratio p of t to R f char- 
acteristic velocity v m ; the terminal densities ??o, 
5 0 and the average densities n, 5 of Fc I atoms. 

The most important result is that using equa- 
tions (16) to (18) , which describe the combination 
of self-absorption and changing volume of radi- 
ating gas in a simple model, we are able to put 
the spectral data on Fe I lines in agreement with 
the luminous efficiency of the total Fe I meteor 
radiation. If temperatures are close to 5500° Iv, 
the values of q are from 0.010 to 0.068, depending 
partly on the total brightness. Extrapolating to 
fainter meteors, there is clear evidence that the 
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effect of a change of the radiating volume with E x 
potentials may be negligible for meteors fainter 
than about — 2 or —1 absolute magnitude. On 
the other hand, the self-absorption effect is 
important even for fainter meteors. If \vc compare 
multiple! 4 with multiplct 15 for the fireballs of 
table 1, and also for the relatively fainter object of 

— 1.5 magnitude published by Harvey (1970), we 
see that the relative self-absorption in bright lines 
of multiplet 4 decreases from about A log I— 1.5 
at -10 magnitude to about A log 7 = 0.8 at 0 
magnitude. Extrapolating, we would expect 
optically thin radiation in multiplet 4 for fainter 
meteors than +5 magnitude, but possibly starting 
with even fainter objects. As multiplet 4 is one of 
the most affected by self-absorption, the meteor 
radiation in the 0 to +5 magnitude interval may 
be optically thin for the majority of Fe I lines, 
that is lines with E^O and sufficiently small gf 
values. But this very preliminary statement 
should be checked in each particular case in 
future studies. 

The resulting relaxation times are in the interval 

— 2.58<log t R < — 1.15, most of the values being 
quite close to an average value of log ^=—1.7. 
This is independent of the actual brightness or 
any other changing parameters. Thus the relaxa- 
tion time of 0.02 s seems to be a typical value for 
30km/s fireballs in the —4 to —12 absolute 


magnitude interval. The light of fireballs is emitted 
during a relatively long relaxation time of the 
gas. This is essential ly the same result as found by 
Harvey (1970) for the spectrum of a —1.5 mag- 
nitude meteor. 

Even if we are working with the absolute values 
of light intensities, it is an advantage that the 
computed values of the relaxation time are not 
dependent on absolute photometry £see equation 
(19) J. They reflect the ratio of the luminous 
efficiency factor to absolute <?/-values, using 
relative photometry for comparison by means of 
the curve of growth. 

The resulting length of the radiating volume is 
between 10 and 70 m, and it clearly depends on 
the actual brightness of the fireball. Extrapolating, 
a 10-m long volume of radiating gas is typical for 
about — 1 magnitude and 100 m is typical for 
about —12 magnitude. It should be mentioned 
that the exceptional spectrum, S 1132, has values 
of t R and t smaller than for the other meteors. 

The diameter 2 R of the radiating volume (for 
Ei = 0) is between 0.4 m and 2.6 m and clearly 
depends on the brightness of the object at the 
point examined (see fig. 6) . A diameter of 28 cm 
seems to be typical for 0 magnitude and a diameter 
of 2 m typical for — 10 absolute magnitude. All 
these values look reasonable. The values of the 
characteristic velocity were computed assuming 




EVIDENCE FROM SPECTRA OF BRIGHT FIREBALLS 


97 



that the atoms moved only along the direction of 
the meteor motion. The resulting characteristic 
velocities, v mj of several km/s in most cases, are 
also what one would expect as the zigzag path of 
the atom means a greater actual velocity than v m . 
The resulting densities of the iron atoms, if 
compared with the standard air densities (19G2) 
at corresponding heights, arc also quite reasonable. 

I tried also to use an exponential change of l 
with the upper potential E 2 . Expressions analogous 
to equations (10), (17) and (18) were derived, 
but in this case the numerical solutions resulted in 
excitation temperatures too high to be admissible 
for Fe I atoms. Very high densities of the radiating 
gas relative to the surrounding air point to a great 
discrepancy in this case. It is clear that the as- 


sumption of a dependence of the radiating volume 
on E'i potentials is not able to explain the dis- 
crepancy between the total light and the spectral 
data. 

Values in table 2 should be considered as one set 
of possible values only. We chose the ratio p as a 
constant value. If p is taken as a function of E\ 
and of time, one can obtain solutions w hich fit the 
observed points even better, but the number of 
free parameters is then too high. Also other models 
than that published in this paper can be used, and 
thus slightly different expressions for the x- and 
y-axes of the curve of growth and for the x-axis of 
the Boltzmann distribution, can be obtained. 
Perhaps some solutions may be found better 
fitting the observational points with such models. 
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Table 2.— Results for Fe I from Emission Curve-of -Growth With Respect to Changes of the Volume of Radiating Gas with 

the Ei potential 


Spectrum 

Point* 

h 

(km) 

b log p 
(g/cm 3 ) 

V 

(km/ 

s) 

Fe I total light 
curve data 

Fe I emission-curvc-of-growth data 
for Ei = 0 

0 4T pan 

Part 
from all 
radiation 
(%) 

d dm/ 
dt 
(g/s) 

No. 

of 

obs ® 

T 

(°K) 

Q 

log N 
(Ei=0) 

S 6 


76 

-7.4 

29 

-9.98 

64 

3.31 

88 

5540 ±160 

0.041 

23.70±0. 13 

S 526 

D 

73.1 

-7.24 

28.3 

-6.41 

53 

1.90 

37 

5500±270 

0.045 

22.78±0. 19 


E 

72.7 

-7.22 

28.1 

-7.42 

37 

2.31 

57 

5510 ±240 

0.051 

22.90±0. 17 


F 

72.1 

-7.18 

28.0 

-7.71 

49 

2.43 

60 

5510±230 

0.064 

22.84 ±0. 17 


G 

71.7 

-7.16 

27.9 

-8.49 

35 

2.74 

79 

5510±200 

0.068 

23. 10±0. 16 


II 

71.2 

-7.13 

27.8 

-7.24 

51 

2.25 

56 

5500 ±220 

0.056 

22.66±0. 16 


I 

70.6 

-7.09 

27.6 

-4.70 

43 

1.24 

16 

5460 ±190 

0.031 

21 .56±0. 13 


J 

70.1 

-7.06 

27.5 

-8.55 

37 

2.78 

76 

5520 ±230 

0.060 

23. 10±0. 17 


K 

69.2 

-7.01 

27.2 

-6.15 

46 

1.83 

38 

5470 ±220 

0.031 

22.49±0. 15 


(M 

68.6 

-6.98 

26.9 

-3.75 

36 

0.88 

10 

5500±240 

0.010 

20.33±0. 16) 

S 912 

D 

77.2 

-7.51 

31.7 

-7.32 

100 

2.17 

22 

5500±450 

0.012 

22.42±0.29 


E 

76.3 

-7.45 

31.6 

-7.76 

100 

2.34 

30 

5520±410 

0.024 

22.70±0.27 


F 

75.7 

-7.41 

31.5 

-8.32 

100 

2.57 

38 

5470 ±370 

0.024 

22. 81 ±0.25 


G 

74.8 

-7.35 

31.3 

-8.68 

100 

2.72 

50 

5470 ±350 

0.034 

22.90±0.24 


H 

74.0 

-7.30 

30.8 

-10.24 

95 

3.36 

115 

5520 ±190 

0.046 

23.67 ±0.15 


I 

73.9 

-7.29 

30.6 

-10.10 

100 

3.31 

108 

5520 ±180 

0.044 

23 . 57 ±0 . 14 


J 

73.5 

-7.27 

30.0 

-10.22 

100 

3.37 

115 

5500±190 

0.046 

23.77±0. 14 

S 1132 


77.0 

-7.49 

32.1 

-12.00 

67 

4.03 

75 

5490 ±130 

0.055 

23.65±0. 12 


But the main purpose of this paper is to demon- 
strate that it is possible to use a simple model of 
fireball radiation which agrees well with both the 
total light efficiency and the spectral distribution 
of the light in absolute values. And this is ade- 
quately verified by the values of table 2. 


RESULTS FOR OTHER ELEMENTS THAN FE I 

With a few exceptions, the whole procedure 
used for Fe I lines could not be repeated for lines 
of other elements. The main reason is the low 
number of lines available for an independent 
solution. 

A solution for 7 Na I lines of spectrum S 1132 
was possible with the result T 7 = 5070° ±170° K 
and Fjj m =1.4 for # = 0.05. The curve of growth 
starts from the region of strong self-absorption 
and continues to the damping region. Thus no 
reliable value of total number of Na I atoms inside 
the corresponding radiation volume could be 


determined. But even if this is not the case, the 
unknown value of the luminous efficiency factor of 
sodium prevents one from estimating the relaxa- 
tion time of Na I. The dimensions of the radiating 
volume of Na I were certainly less than that of 
Fe I as is clear from Y n m being 1.2 less in log for 
Na I than for Fe I. Thus about 4 times less linear 
dimensions of the radiation volume are found for 
Na I atoms than for Fe I. 

A solution for 7 Ca I lines of spectrum S 1132 
was also possible with the result T= 4750° K±300° 
and Y u m = 2.G for g = 0.05. The dimensions of the 
radiating volume are the same as for iron in this 
case. A slightly better position of the emission 
curve of growth than for Na I makes possible an 
estimation of the number of Ca I atoms in the 
whole radiating volume: log N = 20.8 and log m G = 
— 1.4, which is about 0.1 percent of the Fe I 
gas mass. 

A solution for 10 multiplcts of O I of spectrum 
S 1132 was possible with a rough result of tem- 
perature 19 000° ± 4000° K and Fii m = 0.6 for 
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Table 2.— Results for Fe I from Emission Curvc-of-Growth with Respect to Changes of the Volume of Radiating Gas with the 

Ei Potential — Continued 


Spectrum 

Point* 


log rriG 
(g) 

f Turn 
(cm 2 ) 

S 6 


1.67 

2.4 

S 526 

T) 

0.75 

1.2 


E 

0.87 

1.5 


F 

0.81 

1.7 


G 

1.07 

1.9 


II 

0.63 

1.5 


I 

0.47 

1.15 


J 

1.07 

1.8 


K 

0.46 

1.1 


(M 

— 1.70 

1.1 

S 912 

D 

0.39 

2.0 


E 

0.67 

2.1 


F 

0.78 

2.2 


G 

0.87 

2.3 


H 

1.64 

2.5 


I 

1.54 

2.4 


J 

1.74 

2.4 

S 1132 


1.62 

2.6 


Fe I cmission-curve-of -growth data for Fi-0 


log Ir 
(s) 

log p 

log R 
(cm) 

log t 
(cm) 

log Vrn 
(cm/s) 

log no 
(cm -3 ) 

log Jo 
(g/cm 3 ) 

log n 
(cm -3 ) 

log 5 
(g/cm 3 ) 

-1.64 

1.81 

1.95 

3.75 

5.39 

14.82 

-7.21 

15.55 

-6.48 

-1.15 

1.81 

1.73 

3.54 

4.69 

14.56 

-7.47 

15.28 

-6.75 

-1.44 

1.81 

1.73 

3.54 

4.98 

14.67 

-7.36 

15.39 

-0.64 

-1.62 

1.81 

1.74 

3.55 

5.17 

14.59 

-7.44 

15.31 

-6.72 

-1.67 

1.81 

1.79 

3.60 

5.27 

14.70 

-7.33 

15.43 

-6.60 

-1.62 

1.81 

1.68 

3.48 

5.10 

14.60 

-7,43 

15.33 

-6.70 

-1.71 

1.81 

1.53 

3.33 

5.04 

13.94 

-8.09 

14.67 

-7.36 

-1.71 

1.81 

1.74 

3.55 

5.26 

14.84 

-7.19 

15.56 

-6.47 

— 1.37 

1.81 

1.62 

3.43 

4.80 

14.59 

-7.44 

15,31 

-6.72 

-2.58 

1.81 

1.22 

3.03 

5.61 

13.64 1 

-8.39 

14,36 

-7.67) 

-1.78 

1.81 

1.83 

3.63 

5,41 

13.91 

-8.12 

14,64 

-7.39 

-1.67 

1.81 

1.89 

3.70 

5.37 

13.99 

-8.04 

14.72 

-7.31 

-1.79 

1.81 

1.89 

3.69 

5.49 

14.12 

-7.91 

14.84 

-7.19 

-1.85 

1.81 

1.90 

3.71 

5.56 

14.17 

-7.86 

14.89 

-7.14 

-1.72 

1.81 

2.01 

3,82 

5.54 

14.60 

-7.43 

15.33 

-6.70 

-1.77 

1.81 

1.96 

3.77 

5.54 

14.65 

-7.38 

15.38 

— 6.65 

-1.63 

1.81 

2.01 

3.81 

5.45 

14,72 

-7.31 

15.44 

— 6.59 

-2.41 

1.20 

2.12 

3.32 

5.73 

14.86 

-7.17 

15.59 

-6.44 


a Notation of points: see the corresponding paper. 

b infT n niT densities from the TT,S. standard Atmosphere^ 196*.. ... , i 

e d/pal is close to absolute panchromatic magnitude, found by comparing stars with meteor spectrum by n 

of panchromatic emulsion. 

d log Tpan = — 1 1.94. . . . . 

e No. of obs, is the number of Fe I spectral lines used for the solution. 

* Fi im = (log(//(XBx)+4^0Ki) limit for optically thick radiation. 


g = 0.O5. The linear dimensions of the radiating 
volume of O I atoms are 10 times less than that of 
Fe I. The number of O I atoms inside this volume 
is log A r ~20.4, giving log w<?« -2.2 (g) . A rough 
estimate of the corresponding atoms of Fel 
inside the same small volume gives about 1 to 
10 percent of O I mass relative to iron mass. 

No other solutions or even rough estimates for 
all the detailed meteor spectra of table 1 were 
found. In my previous papers, I used a method of 
estimating the abundances of other elements than 
iron from the curve of growth constructed for 
Fe I and applied to a small number of lines, the 
only ones available for the corresponding element. 
But this is possible only on the assumption that 
(a) the excitation temperature T is the same as 
that for Fe I; (b) the dimensions of the radiating 
volume and the characteristic velocity are the 
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Figure 6. — The dependence of the radiating- volume 
radius (log R) of Fe I on the total actual brightness of 
Fe I as derived from the spectral analysis of fireballs 
with 30 km/s velocity. 
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same as that for Fe I; (e) the damping constant is 
the same as that for iron. The independent 
solutions given above even if only rough guesses, 
point clearly to the non-validity of all three 
assumptions. Nevertheless, I tried to find some 
abundances of other elements than iron by as- 
suming that the prepositions (a), (b) and (c) arc 
valid. The results were very frustrating. 

In summary, it does not seem possible to 
determine any realistic abundances of other 
elements than Fe I for fireballs from —4 to — 12 
magnitude. The radiation we receive from these 
fireballs originates mainly from the effective 
surface of the radiating volume, and we do not see 
most of the light originating inside the volume. 
Thus we are able only to compare some total 
numbers of atoms radiating from a completely 
unknown surface of the radiating volume. This 
also strongly depends on what line of what 
element we observe. The situation should be some- 
what better for meteors as faint as 0 magnitude, 
but this is necessary to check by construction of 
the emission curve-of-growth for each particular 
case. As mentioned before, the optical thickness 
still seems to be an important factor for most 
spectral lines of 0 magnitude meteors. 

CONCLUSIONS 

(a) The radiation of fireballs (v^30 km/s, 
70 to 80 km, M pan « — 4 to —12 mag) is 

strongly affected by self-absorption and to some 
extent by damping. 

(b) Self-absorption should still be quite signifi- 
cant at 0 absolute magnitude. Optically thin 
radiation may be expected for meteors fainter 
than +5 absolute magnitude. 

(c) If iron lines of short-duration flares are 
studied by the emission curve of growth and the 
resulting absolute abundances compared with the 
total luminosity by means of the luminosity 
equation factor, the resulting relaxation time is 
longer than the duration of the flare. 

(d) If we assume that the dimensions of the 
volume of radiating gas decrease with increasing 
lower potential E h a relaxation time shorter than 
the duration of the flare can be found. In this way, 
the spectral data on Fe I lines can be made to 
agree with the luminous efficiency of the total 
Fe I meteor radiation. 


(e) Parameter q f relating the dimensions of the 
radiating volume to the lower potential E u was 
found partly dependent on meteor brightness. The 
effect of change of the radiating volume with 
potential E x is negligible for meteors fainter than 
about —2 to — 1 absolute magnitude. 

(f) Excitation temperatures T — SSOCPK were 
found to be an optimum solution for Fe I in all 
the fireball spectra studied. This is effectively a 
lower limit, since somewhat higher values are 
also possible. 

(g) A relaxation time of 0.02 s was found as a 
typical value, independent of the actual brightness 
of the fireball. The light of fireballs is emitted 
during a relatively long relaxation time, which is 
many orders of magnitude longer than the time 
necessary for spontaneous radiation of excited 
Fe I atoms. 

(h) The dimensions of the radiating gas volume 
depend on the fireball brightness. Fe I gas: 
0.3X9 m at 0 absolute magnitude, 2X60 m at 
— 10 absolute magnitude. 

(i) The terminal densities at the end of the 
radiat ing volume of iron vapor roughly correspond 
to the standard air densities at the corresponding 
heights. 

(j) Practically, it was not possible to deter- 
mine any realistic abundances of other elements 
than iron for the fireballs studied. A few lines are 
available for the same independent procedure as 
used for Fe I. The curve of growth derived for 
Fe I cannot be used for other elements. The 
radiation we receive from fireballs originates 
mainly from the effective surface of the radiating 
volume, which prevents a sufficiently accurate 
comparison for various elements of the total 
numbers of atoms inside the whole volume. 

(k) The relaxation time and other problems 
contained in this paper could be studied much 
easier, if meteor spectra with a very short exposure 
time (instantaneous exposure) were available. 
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converting from magnitudes to the absolute scale of I in 
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erg/s. In the phrase “3.59 magnitude = 2.73 times less in- 
tensity than 0 magnitude” the value 2.73 should be 27.3. 

ATI data in my 19G6 paper given in magnitudes or “bolo- 
metric” magnitudes are OK, all data given in log / 
(erg/s) must be corrected by subtraction of 1.00. Table 6 
starts correctly with 11.91 for log I of point H of line No. 1 
and so on till 12.85 for log I of point J of line No. 189. In 
figure 2: the y-axis must be changed from 12,5 to 11.5 
and from 13.5 to 12,5, 

The results in my 1967 paper were derived using the 
wrong numbers mentioned above, but the systematic error 
in the absolute values of the oscillator strengths of Fe I 
acted in the opposite direction. In any case, the problem 
of the analysis of meteor spectra has been studied from 
new and more general aspects in this paper, using the cor- 
rect log I for S 912. 

In my previous papers, I used the notation “bolomctric” 
magnitudes, to indicate that bolometric magnitudes of 
stars were used for computation of the absolute e.g.s. scale 
of intensities. However, the meteor spectral lines were 
compared with stars using a panchromatic emulsion, and 
it is usual to speak about “panchromatic” magnitudes in 
such a case, particularly, if A O comparison stars are used. 

LIST OF MATHEMATICAL SYMBOLS 

a Angle between the line of vision and 

the meteor motion 

B\ Planck's function related to wave- 

length units 

c Velocity of light 

5 0 Terminal density of the radiating vol- 

ume (g/cm 3 ) 

3 Average density at the end of the 

radiating volume (g/ cm 3 ) 
e Charge of electron 

Ei The lower level potential (V) 

E 2 The upper level potential (V) 

F(X) Function giving the curve of growth 

/ Absorption oscillator strength 

g=gi Statistical weight of the lower level 

G 2 Total statistical weight of the upper 

term 

k Planck's constant 

Xi Energy of the lower level 

I Intensity of optical meteor radiation 

I Q Optically thin intensity 

k Luminous efficiency of the radiating 

gas per unit gas mass 
k B Boltzmann's constant 

l Length of the radiating volume 

A Wavelength 


m 
m a 


m e 

m a 

niz 


no 

n 

N 

N( term) 

V 

Q 


R 

t 

tf 

Ir 

T 


Meteor mass 
Initial meteor mass 

Total gas mass, which entered into the 
radiating volume from the be- 
ginning of the meteor luminous 
trajectory till time t 
Mass of electron 

Total gas mass inside the radiating 
volume 

Total gas mass which left the radiating 
volume from the beginning of the 
meteor's luminous trajectory till 
time t 

Mass of one atom 

Terminal number of atoms per cm 3 at 
the end of the radiating volume 
Average number of atoms per cm 3 of 
the radiating volume 
Total number of atoms of one element 
in one state inside the radiating 
volume 

Total number of atoms being in all 
upper levels of one term inside the 
radiating volume 
t/R 

Exponential factor giving the decrease 
of the length of the radiating volume 
with increasing potential E\ } defini- 
tion (14). 

Radius of the radiating cylindrical 
volume 
Time 

Duration of a flare 

Relaxation time of meteor radiation: 
time spent by an average atom 
inside the volume of radiating gas 
Excitation temperature 
Luminous efficiency factor 


5040 



Boltzmann temperature 


U Partition function 

v Meteor velocity 

v m Characteristic velocity 

X x-axis of the curve of growth 

Yum Limit value of y-axis of the curve of 

growth for very strong self-absorp- 
tion 
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